The magnetic field in the NGC 2024 FIR 5 dense core 
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ABSTRACT 



We used the Submillimeter Array (SMA) to observe the thermal polarized 
dust emission from the protostellar source NGC 2024 FIR 5. The polarized 
emission outlines a partial hourglass morphology for the plane-of- sky component 
of the core magnetic field. Our data are consistent with previous BIMA maps, 
and the overall magnetic field geometries obtained with both instruments are 
similar. We resolve the main core into two components, FIR 5A and FIR 5B. A 
possible explanation for the asymmetrical field lies in depolarization effects due 
to the lack of internal heating from FIR 5B source, which may be in a prestellar 
evolutionary state. The field strength was estimated to be 2.2 mG, in agreement 
with previous BIMA data. We discuss the influence of a nearby Hll region over 
the field lines at scales of ~ 0.01 pc. Although the hot component is probably 
compressing the molecular gas where the dust core is embedded, it is unlikely that 
the radiation pressure exceeds the magnetic tension. Finally, a complex outflow 
morphology is observed in CO (3 — > 2) maps. Unlike previous maps, several 
features associated with dust condensations other than FIR 5 are detected. 

Subject headings: ISM: individual (NGC 2024) - ISM: magnetic fields - Polarization 
- Stars: formation - Techniques: interferometric. 
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Introduction 



Understanding the evolution of molecular clouds and protostellar cores is one of the 
outstanding concerns of modern astrophysics. Particularly, efforts are concentrated in 
determining which physical agents are mainly responsible for controlling the dynamical 
properties of the dense cores. It is widely accepted by the astronomical community that 



magn etic fields must 



cores (IShu et al. 



De taken into account in evolutionary models of collapsing protostellar 



19991 ) . Although some the ories claim that turbulent supersonic flows driy es 



star formation in the interstellar medium (jElmegreen fc Scalo 



200 



Mac Low fc Klessen 



2004 ). others demonstrate that the ambipolar diffusion collapse th eory reproduces proper ly 



observed molecular cloud lifet imes and star formation timescales f ITassis fc Mouschovias 



2004 



Mouschovias et al. 



2006|). 



One way to resolve these open issues in star forming theory is to increase the number 
of high-quality observations which resolve the core collapse structure. Sampling several 
protostellar cores with distinct physical properties can provide better constraints to improve 
simulations. Particularly, the number of observations of magnetic fields in molecular clouds 
and dense cores has been increasing rapidly with the advent of new instruments with high 
sensitivity. In polarimetry, i t is globally accepted that non-spherical dust grains are aligned 



perpendicular to field lines ( 



continuum emission ( iHildebrand 



Davis fc Greenstein 



1988 



195ll ) producing linearly polarized thermal 



Which mechanism ma inly contribute s to the 



alignment of interstellar dust grains is still a matter of debate (ILazarian 



very recently 



Hoang fc Lazarianl ( 12008 



20071 ). However, 



20091 ) have successfully modeled the polarization 



by radiative torques propelled by anisotropic radiation fluxes. Those torques act to align 
spinning non-spherical dust grains with their largest moment of inertia axis parallel to the 
field lines. The polarized flux is usually only a small fraction of the total intensity (usually 
a few percent) and for this reason the study of the magnetic field is highly limited by the 
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instrumental sensitivity. 



Cold dust emits mainly at far-IR and submm wavelengths. In the submm regime, 
the emission is optically thin and, therefore, it is not affected by scattering or absorption. 
For this reason, the SMA has been extensively used to study thermal emission from dust 
and cool gas. Several authors reported polarization observations of different classes of 
protostella r cores. The textb ook case is the low mass young stellar system NGC 1333 



IRAS 4A flGirart et al 



20061 ). The supercritical state is reflected in the SMA polarization 



maps which indicates a clear hourglass morphology for the plane-of-sky magnetic field 



component in a physical scale of 300-1000 AU. This re markable result not on 



predicted by theo ries of collapse of magnetized clouds ([Fiedler &: Mouschovias 



y was 



1993 



Galli k Shu 



(IShu et al 



2006 



1993 ) but is commonly used to test mode 



Goncalves et al 



2008 



Rao et al. 



s of low-mass collapsing cores 



20091 1 . In the regime of high mass 



protostars, recent investigations performed with the SMA also provided observational 
constraints on the physics involved during the core collapse stage. Two recent works 
exemplify distinct magnetic field features within this class of objects. The polarimetric 
properties of G5.89-0.3 9 are consistent w ith a complex, less ordered field likely disturbed 



by an ionization front (ITang et al. 



20091 1. while the hourglass morphology expected for 
magnetically supported regi mes was observed a t large physical scales (~ 10 4 AU) for the 



protostellar core G31. 41+03 (IGirart et al. 



20091 ). Despite the distinct energy balance and 



timescales of the two regimes (low and high mass), both results imply that the magnetic 
support must not be ignored in the models. 

NGC 2024 is the most active star forming region in the Orion B giant molecular cloud. 
The gas structure in this region has an ionized component surrounded by a background 
dense molecular ridge and a foreground dust and molecular component visible in the optical 
images as a dark lane. Recently, the new ESO telescope VISTA (Visible and Infrared 
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Survey Telescope for Astronomy) released a high sensitivity near-infrared image of NGC 
2024 (Figured]). In this large scale view, the foreground dust lane which optically obscures 
the Hll region is almost transparent. Scattered light produced by the ionization front is 
seen as bright emissi on in the top of the image, and a cluster of hot young stars is revealed. 



Kandori et al 



(120071 ) used near-infrared polarimetry to study the scattered light from the 
Hll region. Several reflection nebulae associated with young stellar objects (YSO) were 
discovered in their polarization maps. The overall centro-symmetric pattern suggests that 
the ionizing source is IRS 2b, a massive star located 5" north-w est of IRS 2 , in a greement 
with a previous near-infrared photometric study carried out by 



Bik et al 



fl2003|). The 



su bmillimeter continuum em ission arising from the dense molecular ridge was first observed 



by 



Mezger et al. 



(11988 



1992j). Several far infrared cores (so the acronym "FIR") at distinct 
evolutionary stages were identified and catalogued in a North-South (NS) distribution. In 
fact, the chain of FIR cores could have been genera ted by the interaction betw een the nearby 
Hll region and the surrounding molecular cloud. iFukuda fc Hanawal (120001 ) performed 
numerical simulations of sequential star formation trigged by an expanding Hll region near 
a filamentary cloud. In their models, isothermal expansion and magnetohydrodynamic 
effects are considered. Their simulations preview that a chain of cores is formed from this 
interaction, each pair of cores belonging to a distinct generation though. Comparison 
between model and the dynamical parameters observed in NGC 2024 are quite good. In 
particular, they state that FIR 4 and FIR 5 belong to the first generation of cores, what is 
confirmed by the observed dynamical ages of their outflows. In this paper, we center our 
discussion on FIR 5, the brightest an d most evolved of them, with a strong and collimated 



unipolar outflow 



Wiesemeyer et al. 



Richer et al. 



19921 ). Continuum observations at 3 mm performed by 



19971 ) suggest that FIR 5 is a do uble core embedd ed in an envelope. 



However, higher angular resolution observations from 



Lai et al 



((2002) (LCGR02 from now 



on) resolved the dust emission in one strong component surrounded by several weaker 
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components in a radius of a few arcseconds. 



Several aut 



Crutcher et al. 



lors have conducted polarimetric investigations toward FIR 5. 



( 119991 ) used the Very Large Array (VLA) to carry out Zeeman ob- 
servations of OH and Hi absorption lines in order to trace the line-of-sight (LOS) component 
of the magnetic field. These authors found a LOS field gradient of ~ 100 fiG across the 
northeast-southwest di rection. Dust emi s sion p olari zation in the surroun dings of FIR 5 was 



mapped in 100 /mi by 



Hildebrand et al 



fll995f ) and 



Dotson et al 



Matthews et al. 



(120001) with the Kuiper 



((20021) used the SCUBA 



Airborne Observatory. At longer wavelengths, 
polarimeter to observe the 850 /im emission with at the James Clerk Maxwell Telescope 
(JCMT) and obtained polarization patterns consistent with those derived with the 100 fim 
data. Their single-dish dust polarization maps trace a relatively ordered magnetic field 
along the ridge of emission containing FIR 4 and FIR 5 (Fig. 4 in their paper). Based on 
the spatial coverage of their observations, this ordered field must extend o ver a distance of 



at leas t ~ 0.5 pc (for the assumed distance of 415 pc to the Orion B cloud). 



Matthews et al. 



( 120021 ) modeled this field using a helical-field geometry threading a curved filament, since 



this configuration suited fairly to a 2-dimensional projection accordingly to the SCUBA 



maps. However, they 



data of 



Crutcher et al 



bund that this geometry is not consistent with the LOS Zeeman 



(119991 ) because no reversed fields are seen at both sides of the chain 
of cores. Instead, those authors offer another interpretation based on a compression zone 
created due to the expansion of the foreground ionization front. In this scenario, the 
magnetic field lines are stretched around the ridge of dense cores i n a physical mo r pholo gy 



consistent with LOS field gradient observed in the Zeeman data of 



Crutcher et al. 



(119991). 



Concerning the local field associated with FIR 5, the work of Lai et al. (2002, hereafter 
LCGR02) has the best resolution for the dust continuum emission so far. These authors 
used the Berkeley-Illinois-Maryland- Association (BIMA array) interferometer and obtained 
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an angular resolution of 2'/4 x V'A. The polarized flux of the BIMA maps extends in a N-S 
direction, perpendicular to the putative protostellar disk. The corresponding field lines 
were fitted with a geometric model consisting of a set of concentric parabolas, indicating 
that the polarized flux trace a partial hourglass morphology for the magnetic field. In this 
paper, we report SMA dust continuum polarization toward FIR 5. The higher sensitivity 
of this instrument provides new information on the detailed field morphology of the FIR 5 
core. 



Observations and Data Reduction 



The high angular resolution of the SMA allows us to trace the thermal emission of 
dust grains at physical scales of few hundred astronomical unitj* (for objects in the Orion 
molecular cloud complex) and, therefore, is able to spatia lly resolve compact dust cores. 



A detailed description of SMA is given in 



Ho et al. 



( 120041 ) . The observations were carried 



out in 2007 November 24 and December 19 with the SMA in its compact configuration. 
The number of antennas available for the observations were 7 and 6, respectively. The 
atmospheric opacity at 225 GHz was 0.11 and 0.07 for the first and second day, respectively 
(values measured by the Caltech Submillimeter Observatory tau meter). Observations were 
done in the 345 GHz atmospheric window, what corresponds to a wavelength of 870 jum. 
The SMA receivers operate in two sidebands separated by ~ 10 GHz. The central observed 
frequencies for the lower and upper side bands were 336.5 GHz and 346.5 GHz, respectively. 
The SMA correlator had a bandwidth of 1.9 GHz (for each sideband) divided in 24 "chunks" 



x The Submillimeter Array is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded 
by the Smithsonian Institution and the Academia Sinica. 
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of 128 channels each. In total, the full-band spectrum contains 3072 channels for each 
sideband and a spectral resolution of 0.62 MHz, which corresponds to a velocity resolution 
of 0.7 km s _1 . SMA receivers are single linearly polarized. By using a quarter-wave plate in 
front of each receiver, the incoming radiation is converted into circular polarization (L, R). 
The SMA correlator combines the signal into circular polarization vectors: RR, LL, RL, 
LR. In order to obtain the full four Stokes parameters for all the baselines, the visibilities 
have to be averaged on a time scale of 5 minutes. A description of the SMA polarimeter 
an d the discussio n of th e me thodology (both hardw are and software aspects) are available 



in 



Marrone et al. 



O2006|) and 



Marrone fc Raol ( 120081 ) 



The phase center (a 2 ooo = 05 h 41 m 44!3, <5 2 ooo = -01°55'4(X / 8) was set according to the 
peak of emission obtained for FIR 5 in LCGR02. Uranus and Titan were observed as flux 
calibrators in both tracks. The resulting visibility function for each calibrator is consistent 
with the expected flux estimated by the SMA Planetary Visibility Function Calculator 
during the observing runs. The quasar J0528+134 was used as the gain calibrator. The 
quasar 3c454.3 was used for bandpass and polarization calibration. The first track had a 
much better parallactic angle coverage for 3c454.3 than the second track, thus the 3c454.3 
data from the first track were used to solve for the instrumental polarization or "leakages" . 
The minimum and maximum UV distance for both tracks was 16 and 88 kA, respectively. 
Antenna 3 was used only in the second track, so no leakage solution could be derived. 
Thus, antenna 3 was not used to obtain Stokes Q and U maps. After the calibration steps, 
data from upper and lower sidebands for each track were synthesized into a single data 
set. Calibrated visibilities for each track were combined into a final data set. Removal of 
continuum contamination from the line data set was done. The main contribution arose 
from the CO (3 — > 2) transition at the chunk ^ 4 of the upper sideband (~ 345.76 GHzfl 



2 Since our interest in the line data set concerns only Stokes I emission, antenna 3 is 
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All the calibration and reduction steps were done with MIRIAD configured for SMA 



data ( jWright Sz Saultlll993l ). The science target was strong enough and self-calibration was 
performed in order to increase the signal-to-noise ratio in the final maps. Imaging of the 
Stokes parameters I, Q and U was performed. Maps of polarized intensity (Ip), polarized 
fraction (P) and position of polarization angles (6) were obtained by combining Q and U 
images in such a way that P = = ^® an d 9 = |tan _1 (^). The resulting synthesized 
beam of Stokes I maps has 2'.'45 x 1'.'48, with a position angle (PA) of —39.8°. Tabled] 
summarizes the technical parameters of continuum and line observations. 



Results 



3.1. Dust Continuum Emission 



Figure [2] shows the contour map of the 878 /im dust emission in FIR 5 obtained with 
a quasi-uniform weighting (a robust of —1), which provides a better angular resolution of 
l'/96 x 1'.'41. The overall submillimeter emission resembles the 1.3 mm dust continuum maps 
obtained with BIMA by LCGR02, although the latter has a rms a factor of 3 lower. Our 

allow us to only measure structures 



Palau et al 



fboinh ). In LCGR02, they 



observations (with shortest baseline equal to 16 kA~ 
smaller than ~ 5.7 arcseconds (see Equation A. 5 of [ 
find extended emission at scales of 5-7 arcseconds with BIMA. The brightest emission arises 
around FIR 5: Main (following LCGR02 notation), which is resolved into two components, 
5 A and 5B. Those comp onents correspond to the double source detected in 3 mm by 



Wiesemeyer et al. 



(119971 ) and indentified as FIR 5-w and FIR 5-e. However, not all the 



fainter sources observed in the FIR 5: Main core of LCGR02 have been detected with the 
SMA. FIR 5: Main appears more extended in the BIMA maps, attributable to the better 



unflagged in the deconvolved CO (3 — > 2) maps. 
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sampling of shorter baselines with the BIMA array. In particular, the N-S direction contains 
emission of the dust condensations LCGR2, LCGR3 and LCGR5 (according to LCGR02 
nomenclature). These sources are missing in our SMA maps probably due to the absence of 
antenna 3 in the deconvolved maps (see section [2]). Antennas 3 and 6 cover a short baseline 
in the UV plane which is parallel to the U axis and close to V = kA. In equatorial 
coordinates, it corresponds to features parallel to the declination axis, and close to the 
phase center. Therefore, by flagging antenna 3 we lost this flux component which should be 
produced by the missing sources. In addition, different visibility sampling between the two 
aperture synthesis telescopes (BIMA sampled shorter baselines) must also be considered. 
We also detect several faint peaks (at the 4 and 7-er level; 1 o = 18 mJy beam -1 ) also seen 
by LCGR02. 

Tables [2] and [3] give the dust emission properties for the two condensations associated 
with FIR 5: Main and for the fainter dust condensations, respectively. The intensity peak 
and the position of the sources were derived using the Miriad task "maxfit" . For the two 
bright sources associated with FIR 5: Main, a two Gaussian fit (using the AIPS's "imfit" 
task) was used to estimate the flux density of each component and its size. The two sources 
appear resolved but in different directions. Thus, source 5A has a full width half maximum 
(FWHM) size of 2"8 x 2"4 elongated close to the NE-SW axis (PA = 67°), whereas source 
5B has a deconvolved FWHM size of 3'.'6 x 2'.'8 but is elongated along the SE-NW direction 
(PA = 124°). 

In order to estimate the column density and mass of the cores, we need to assume a 



1993 



Mangum et al 



1999 



Watanabe fc Mitchell 



2008 



> between 40 and 85 K ( 


H 


o et al. 


Emprechtinger et al. 


2009 


). Here 



adopt a temperature of 60 K. We assume a dust opacity at 878 fim of 1.5 cm 2 g , which 
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approximatel y is the expected value for d ust grains with thin dust mantles at densities of 



10 6 cm 3 (jOssenkopf fc Henning 



19941 ). Using the previous FWHM sizes derived from 
the Gaussian fit, a beam-averaged column density of ~ 4.7 x 10 23 and 2.2 x 10 23 cm -2 for 
sources 5A and 5B were derived, respectively. Similarly, masses for these two components 
are 1.09 and 0.38 M , res pectively. The total m a ss of FIR 5:Main, 1.5 M , is consistent 



with the value derived by 



Chandler &: Carlstroml (119961 ) . 



3.2. Distribution of the polarized flux 

For better sensitivity to the weak polarized emission, maps of Stokes I, Q and U 
were obtained with a robust weight of 0.50. Figure |3] shows the Stokes Q and U emission, 
which have different distributions. The Stokes Q emission arises from a negative compact 
spot about one arcsecond north of source A. The Stokes U is quite extended along FIR5: 
Main, with the brightest emission around source 5A. Source 5B has only weak polarized 
emission at 3-<j level. It is noteworthy that significant positive Stokes U emission appears 
west of source 5A without dust emission associated. However, this spot of polarized 
emission coincides with the BIMA continuum source FIR 5: LCGR 1 (catalogue of 
LCGR02). The non-detection by the SMA could occur because dust emission has been 
resolved out by the interferometer (approximately 30% of the flux is filtered out, see 
section IBTTj) . Thus, we tentatively associated this polarized spot to this source. A cutoff 
of 2-cr (1- a ~ 5 mJy beam -1 ) in polarized intensity (a/Q 2 + U 2 ) is used to obtain the 
linear polarization emission and to derive the position angle in the plane of the sky of the 
polarization vectors. 

The polarization intensity and the polarization fraction in our maps achieve values as 
high as 54 ± 6 mJy beam -1 and 15% ± 2%, respectively, at the northern portions of the 
core, where the polarized emission is brighter. Figure H] shows the dust continuum emission 
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from the protostar overlaid with the dust polarization vectors. Using the position of the 
continuum peak as reference, three main components can be distinguished: a northern 
component, where the highest polarization degrees are obtained, a southwestern component 
and an eastern component offset by ~ 5" from the continuum peak. This distribution is well 
represented in the histogram of polarization angles shown in Figure [H There is a change of 
roughly 90° in the position angles of vectors associated with FIR 5A and the eastern vectors 
associated with FIR 5B. Concerning only vectors associated with FIR 5A, position angles 
have a gradual rotation of approximately 40° from north to south. Table H] summarizes 
our polarization data. Note that the distribution of the polarized flux of the SMA data 
is remarkably consistent with the BIMA maps of LCGR02. Although the structure of 
emission in both the BIMA and SMA data sets has the same ove rall pattern, the l atter has 



Matthews et al. 



fl2002h 



a larger area of polarized flux. Compared to the JCMT maps of 
the mean direction of our SMA polarization field is consistent with the lower resolution 
single-dish data, which do not resolve the structure of FIR 5 and traces a larger physical 
scale associated with the diffuse gas found at the core envelope. 



3.3. CO (3 2) emission 

Our SMA CO (3 — > 2) maps reveal a very complex morphology possibly related to 
multiple outflows. Figure [6] shows the channel maps of the CO (3 —> 2) emission with a 
velocity resolution of ~ 2.1 km s -1 . At blueshifted velocities the emission arises from an 
elongated but wiggling structure in the East- West direction. This blueshifted component 
appears to be associated with FIR 6. The distribution of the molecular gas at the cloud 
systemic velocity (i>lsr — 10 km s _1 ) is basically associated with the FIR 5 main core. 
At redshifted velocities (i'lsr ^ 14 km s _1 ) there are two main elongated features almost 
parallel extending in the North-South direction and observed over a wide range of velocities 
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(up to f lsr — 30 km s 1 ). One of these featu r es is associated wi t 



i the well-know outflow 



powe red by FIR 5A (1 Sanders fc Willner 



1985 



Richer et all 



1992 



Chandler fc Carlstrom 



19961 ) and the other one is located about 10" to the west and seems to arise from FIR 5-sw. 
These two lobes have their brightest emission located near their associated dust components 
(FIR 5A and sw). The emission presents a clumpy morphology, with an average angular 
size of ~ 5", corresponding to a physical size of 0.01 parsecs. It is worth noting that the 
three possible CO high velocity features have no counterpart at the opposite flow velocities. 
Thus, the North-South redshifted lobes have no blueshifted counterpart, and the East-West 
blueshifted lobe does not have a redshifted counterpart. 

Figure [7] shows the Position- Velocity (PV) diagram centered in FIR 5A with a PA 
= 0°.9 (along the brightest redshifted lobe). The outflow powered by FIR 5A has a wide 
distribution of velocities which prevails until ~ 30 km s _1 . An extended spatial distribution 
is observed up to ~ 35" south of the source, although only low velocity components are 
observed at such distances. No blue lobe is seen and only residual emission is measured in 
the northern counterpart. The blue component observed at the offset position of —14" is 
part of the outflow associated with FIR 6. In section I4~4l we provide a detailed discussion 
about the molecular distribution in this region. 



4. Discussion 

4.1. Polarization properties 

In this section, we focus our analysis on the northern and southwestern polarization 
features, which are the brightest components and scientifically more interesting since a less 
uniform pattern is observed. At 2-cr level, these two regions are connected and surround the 
peak of total intensity. From Figure HJ it can be noted that the peak of polarized and total 
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intensities are approximately 2'.'6 apart. Figure [8] shows the dependence of the polarization 
fraction with Stokes I flux and with respect to the distance to source A. In both cases, 
there is a clear depolarization toward the center, where the highest density portions of the 
core are located. The left panel of Figure suggests that the distribution of polarization 
with respect to the Stokes I emission seems to be composed by two subsets: the highest 
polarization fraction data that has a slower growing curve and corresponds to the northern 
component, and the subset with a linear dependence, which arises from the southwestern 
component. The right panel of Figure |8] was produced by performing averaging over 
polarization data for concentric annuli of 0'/4 each. 

The depolarization effect is observed not only at the brightest component, source A, 
but also for the second dust component, source B, represented by a "hole" at r ~ 4'/ 5 in the 
right panel of Figure [SJ Those diagrams are consistent with Figure IU where the polarization 
fraction increases with distance from the peak of emission, but there is a lack of overall 
polarized emission toward source B. 



The depolarization observed at high er values of Sto 



effect observed at different wavelengths ( iGoodman et ai- 



res I seems to be part of a g lobal 



1995 



Lazarian et al. 



1993). In 



the mm/submm range, this phenomenon was also observed in the B IMA data publis 



red 



by LCGR02, as well as far- infrared observations with single dishes ( ISchleuning 



1998 



Matthews et al 



200lf ) . The anti-correlation between Stokes I and polarization fraction 



can be caused by different mechanisms. On one hand, it may be the result of changes 
in the grain structure at higher densities. Those changes may be responsible for a 
decrease in the efficiency of dust grain alignment with respect to the local magnetic field 



( iLazarian fc Hoang 



20071 ). In the case of FIR 5B, the embedded source may be in a very 
early stage of formation, prior even to collapse (since no clear evidence of star-forming 
signatures like outflows has been assigned to it). In this case, the lack of internal infrared 
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radiation could provide no radiative torque to the dust grains and, therefore, no polarized 
flux is observed. Another explanation for this effect could be a twisted magnetic field or 



the superposition of dist inct field directions a 



net polarization degree ([Matthews et al. 



ong the LOS resulting in a reduction of the 



20011 ). Observations at higher angular resolution 



would be necessary to resolve the small scale structure. 



4.2. Magnetic field properties 



In section [T] we briefly introduce the physical mechanisms associated with the alignment 
of dust grains with respect to the magnetic field lines. Although some works propose that 



grain alignm ent co u 



particle flux, 



d be independent of magnetic fields (e. g.: mechanical alignment by 



Goldl ( 119521 )). it has not been proven yet observationally. Dust grains are 



believed to have at least a small fraction of atoms containing magnetic momentum in their 
compositions, so some interaction with the ambient magnetic field is expected. The most 
stable energy state is achieved when the grain longest axis rotates perpendicularly to the 
field lines. Consequently, dust emission polarization vectors as observed in submillimeter 
polarimetry have to be rotated by 90° in order to be parallel to the plane-of-sky (POS) 
component of the magnetic field. The LOS field component adds no information to the 2D 
polarization map because the spinning dust grains produce zero polarization flux. If a strong 
LOS component is expected, a decrease in net polarization flux is observ ed, and alternative 
techni q ues must be used to m easure it (e.g., Zeeman effect observations: 



(|1982j) 



Crutcher et al. 



Troland fc Heiles 



( 119931 )). Therefore, the polarization map of Figure HJ when rotated 
by 90°, traces the projection of the 3D magnetic field morphology on the plane-of-sky (see 
Figure [S]). For FIR 5A, the field geometry is described by curved lines centered on the 
protostellar core. Toward the elongated emission associated with FIR 5B, the field lines are 
parallel to the core's major axis, implying a 90° change in the direction with respect to the 
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FIR 5A mean direction. By relaxing the signal-to-noise level down to 1-er, one can see that 
this change in the magnetic field direction is not abrupt, and an hourglass morphology can 
be roughly derived for the main component (Figure [HI upper right box). Several theoretical 
works have performed 3D simulations of collapsing magnetized clouds. They all agree that 



the PO S projection o 



shape (lOstriker et a 



Girart et al. 



fl2006h 



the magnetic field morpholo gy in those class of objects is a hourglass 



2001 



Goncalves et al. 



Rao et al 



20081 ) . Our results, and many others (e.g. 



(120091 )) provide observational support to these models. 



So far, the CF relation developed by 



Chandrasekhar fc Fermil (119531 ) is still the most 



straightforward method to estimate the plane-of-sky component of the magnetic field. 
Assuming energy equipartition between kinetic and perturbed magnetic energies as 

\pWlos ~ ^5B\ (1) 

(where SVlos is the observational rms velocity along the line-of-sight and p is the average 
density), this method compares the fraction of uniform to random components of the field 
under effects of Alfvenic perturbations (Sv oc 5By/~p) taking into account isotropic velocity 
dispersions. The CF formula uses the dispersion of position of polarization angles and 
molecular line widths as observational inputs for the gas motions in the core. However, 



recent works showed that 



resolutions ( Heitsch et al. 



;his approximation overest imates the magnetic field for coarser 



2001 



Ostriker et al. 



20011 ) . These authors constrained the 



application of this method to data sets with relatively low dispersion of position angles 
(A8 < 25°), which means strong-field cases. By statistical studies of magnetic turbulent 
clouds, these authors showed that the CF formula is accurate only if this condition is 
applied. Using the small angle approximation 5(p ~ SB/ B un if orm , the CF formula can be 
stated as follows: 

B un iform = ^\f^P TT~ i ( 2 ) 



where 5(f) is the angle dispersion. The correction factor £ (~ 0.5) arises from the previously 
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mentioned strong field conditions to which this case applies ( lOstriker et al.ll200lf ) 



Unlike LCGR02, we opted for not applying any geometric model to the observed 
field due to the low statistics of our data set. The observed dispersion in our data (main 
component in the histogram of Figure |5]) is 12.2°. According to 5<f) b s = {^4>int + a e) 1 ^ 2 ^ 
the observed dispersion depends on the intrinsic dispersion S<pi n t plus the measurement 
uncertainty of the position of the polarization angles ag, resulting from the contributions 
of both effects. Since that no geometric models were used to remove the systematic field 
structure, changes on the large-scale field directions are included in the intrinsic dispersion, 
together with turbulent fields due to Alfvenic motions. In our data set, the position angle 
uncertainties average to 7.52°, which gives us an intrinsic dispersion of 9.61°. Some extra 
observational parameters are needed to compute the magnetic field strength with equation 
EJ The average density and therms velocity in FIR 5 can be obtained from previous 
works. lEmprechtinger et all ( 120091 ) modeled the morphology of NGC 2024 based on APEX 
observations of CO isotopologues. The various line profiles obtained for different transitions 
are consistent with a complex structure composed by a Photo Dominated Region (PDR) 
foreground to the molecular gas where the far infrared cores are found. In their models, 
the dense molecular cloud must be warm (75 K) and dense (9 x 10 5 cm" 3 ) to reproduce 
the velocity gradient s observed for d i stinct cloud components. These results agree with 



the previous work of 



Mangum et al. 



fjl999l ). based on formaldehyde observations. These 



authors derived a kinetic temperature of T# > 40 K for FIR 3-7 and estimate densities at 
the same order of magnitude (nu 2 ~ 2 x 10 6 cm -3 ). We adopt wh 2 = 1.5 x 10 6 cm -3 as 
an average value for the density. For the v elocity dispersion, we adopt SVlos of 0.87 ± 



0.03 km s x , which is the value derived by 



Mangum et al. 



( 119991 1 from the formaldehyde 



observations. This m olecule is a good tracer of dense gas, and for the single-dish data of 



Mangum et al 



(Il999l ). it traces the gas kinetic temperature in a scale of ~ 8000 AU, hence 



it is well correlated to the turbulent motions of the core. Finally, applying those inputs 
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to the equation [2], together with the 5<f>i n t previously obtained, we estimate that the POS 
magnetic field strength is 2.2 mG, which is in good agreement with the value estimated in 
LCGR02. The uncertainty in the magnetic field strength is determined mainly by the error 
in the volume density n, which is a factor of ~ 2 due to the distinct assumptions on the 
cloud temperature. This factor implies an uncertainty of 40% for the derived field strength. 
As mentioned earlier, the dispersion used as input in equation [2] carries the combined 
effects of changes on the large-scale field directions plus turbulent motions. In this case, 
the derived field strength is only a lower limit since the angle dispersion is not generated 
purely by Alfvenic motions. On the other hand, beam averaging and line-of-sight effects 
due to field twisting of multiple gas components usually underestimates the real value of 
the turbulent component, and the estimated field strength in this case would be an upper 
limit. So, we can assume that both effects cancel out and 2.2 mG is a fair estimation for 
the POS field strength. 

The mass-to-flux ratio gives information on whether the magnetic field can support 
the cloud against the gravitational collapse and, therefore, it provides clues about the 
evolutionary state of the source. Specifically, this quantity compares the pressure produced 
by an amount of mass M in a magnetic tube of flux $. A critical value, reached when 



the magnet ic pressure is no longer ab 





Nakano & Nakamura 


1978) 


(Crutcher et al. 


1999 


»: 



e to support the gravitational pulling, is given by 



19781 ). Observationally, this parameter is defined by 



A = { ^^ observed = (mN(H 2 )i/5i) x (2tcVG) = 7.6 x 1(T 21 ^^, (3) 

(M/<P) critical B 

where [Mj §) C riticai is the mass-to-flux ratio of an uniform disk where gravity is balanced by 
magnetic pressure, m = 2.8m# allowing for He, A is the cloud area covered by observations, 
N(H 2 ) is the column density and B is the magnetic field strength. Applying the POS 
magnetic field strength obtained in the previous paragraph, B = 2.2 mG, and the column 
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densities derived in section [37T| we estimate a mass-to-flux ratio for FIR 5A of 1.6 (for 
T = 60 K), which corresponds to a core in a supercritical stage. In any case, those 
calculations are restricted to the dust envelope, without taking into account the mass 
contribution of the embedded protostar. We consider that the derived mass-to-flux values 
are only a lower limit for this quantity and therefore it is in agreement with the observed 
star-forming signatures. 

We can also derive the ratio between turbulent and magnetic energies. From the 
autocorrelation function of the polarization position angles, it is possible to measure how 
the dispersion of PA's varies with respect to the distinct length scales within the cloud. 
This function provides an indirect calculation of the turbulent to magnetic energy as 



For the angular dispersion obtained from our sample, 5(J)i nt = 9.61, we compute the turbulent 
to magnetic energy ratio as (3 tU rb = 0.33. This value agrees with the ratio estimated in 
LCGR02, which reinforces that the turbulent motions are magnetically dominated. The 
turbulent-to-magnetic energy ratio found for FIR 5 is consistent with what was measured 
for other low-mass protostellar cor es like NGC 1333 IRAS 4A and IRAS 16293 - 2422 



(Girart et al. 


2006; 


Rao et al. 


2009) 



4.3. Magnetic field around FIR 5A: gravitational pulling or Hll compression? 

In this section, we try to elucidate which mechanism is mainly responsible for 
the observed curved magnetic field morphology in FIR 5. One possibility is that the 
gravitational pulling overcomes the local magnetic support and drags the ionized material 
toward the center, warping the field lines in such a way that they assume an hourglass 
morphology. This is consistent with the previous result that the protostellar core is in the 



flHildebrand et all 120091 ): 




(4) 
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supercritical regime. However, if this is the case then only the hourglass component west 
of FIR 5A is observed. The lack of detected vectors east of FIR 5A could be due to the 
overlap in the line of sight of the dust polarization associated with both FIR 5A and FIR 
5B cores. The magnetic field direction associated with FIR 5B is perpendicular to the FIR 
5A main direction. Alternatively if the two cores are connected, then it could be due to an 
abrupt change in the magnetic field direction. In both cases, the net polarization flux is 
expected to decrease significantly. Another possibility is that the grain alignment efficiency 
associated with source B is smaller. Indeed, Figure [8b (right panel) show that the second 
polarization "hole" matches quite well to the position of FIR 5B. Of course, the cause could 
also be a combination of these possibilities. 

If the tension generated by the magnetic field curvature is produced by the gravitational 
collapse, then we can make a rough estimation of the mass required to produce the observed 
curvature. This magnetic force can be expressed as B 2 / R, where R i s the r adius of curvature 



of the field lines. According to the equations derived by 



Schleuningl ( 119981 ). we have 



M 



100M P 



0J 



B 



lmG 



D 



O.lpc 



R 



0.5pc 



n(H 2 



10 5 cm 



(5) 



where D is the distance of the field lines from the protostar. At D — 1"9 (789 AU) the field 
lines have a radius of curvature R of 17" (7055 AU). At the selected radius of curvature, 
the estimation of magnetic tension force is ~ 10~ 23 dyne cm -3 . Applying these values to 
equation |5l we find that the mass inside the radius of l'/9 is ~ 2.3 M . Although this 
value is almost a factor of two higher than the mass estimation done for FIR 5A in section 



I3.1[ it is within the same order of magnitude of the first estimation, even with the large 
uncertainties in the assumptions of D and the radius of curvature R. This method is an 
alternative approximation to test if gravitational pulling is the responsible for the magnetic 
field curvature. 



Given the situation of the FIR 5 core, the external agents may also interfere in 
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the protostellar physical environment. Previous observations proved that the molecular 



(Barnes et al. 


1989: 


Schulz et al. 


1991: 


Chandler & Carlstrom 


1996: 


Crutc 


ler et al. 


1999 




The distribution of molecular and ionized gas proposed by 


Matthews et al. 


( 


2002 


) for NGC 



2024 (Figure 8 in their paper) has the western portion of the ionization front expanding 
toward the background molecular ridge and stretching the magnetic field lines around the 
ridge of dense cores. At large scales (~ 0.5 pc), this morpholo gy is corroborated no t only 



1999|) but 



by the LOS field obtained from the CN Zeeman observations (jCrutcher et al. 
also by the POS field from the single-dish dust polarization data. At smaller scales (~ 0.02 
pc), this could have an effect of compressing the magnetic field lines, bending them toward 
the east, as observed around the FIR 5 core. In order to check if the radiation pressure 
can be large enough to compress the molecular gas, we have studied the distribution of 
the ionized gas in NGC 2024. For this purpose, we accessed the NRAO Data Archive 
System to search for centimeter emission that could reproduce this morphology. We found 
an extended emission in 6 cm related to the Hll region produced by the star IRS 2b. 
Figure [TU] shows that the hot gas has an extended component to the west and is roughly 
flattened to the south (although slightly curved to the southwest). FIR 5A and FIR 5B, 
indicated as crosses in Figure [TQl lie in the border of the Hll region. The bright southern 
pattern near FIR 5 could trace the compressed ionized gas resulting from the interaction 
between hot /diffuse and cold/denser components. The radiation pressure produced by the 
illuminating star (P ra d) can be calculated by A, where L is the luminosi ty of the ionizati on 



Bik et al. 



(120031 ) 



source, c is the speed of light and A is the area of the expanding shell, 
found that the spectral type of IRS 2b is in the range 08 V-B2 V, which is consistent 
with the intensities of radio continuum and recombination lines observed in the Hll region 



(IKruegel et al 



1982 



Barnes et al 



19891 ). Therefore, we can assume L = 10 L Q , which is 



representative of such a spectral type. A first estimation for the radius of the Hll region 
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was done by ISchraml fc Mezgerl (jl969l ) through low resolution (~ 2') radio observations 
of NGC 2024. These authors measured a radius of 0.2 pc (~ 41 x 10 3 AU) inferred from 
the observed emission. However, from Figure [TD1 the radius of the centimeter emission 
can be fairly estimated in ~ 1', which is approximately the distance between FIR 5 and 
IRS 2b. As a result, the radiation pressure P ra d is calculated as 1.16 xl0~ 8 dyne cm -2 . 
The ionization pressure (N e x T e x k) also accounts for the energy produced by the PDR. 
We assume an electron density of 5.94 xlO 3 cm -3 as derived fro m the emission paramet ers 



Reifensteinetal 



dl970h 



of the centimeter VLA map. The recombination line studies of 
provide an electron temperature of 7200 K for this Hll region. To be conservative, we adopt 
a range of 7200-15000 K for T e . With these parameters, the ionization pressure is estimated 
to range between 5.9 x 10~ 9 and 1.2 x 10~ 8 dyne cm" 2 . On the other hand, the magnetic 
pressure is defined by P mag = J^-, where B is the total magnetic field strength. Since our 
SMA maps provide a two-dimensional picture of the total field, we are able to calculate only 
a lower limit for the magnetic pressure. Therefore, applying the previous equation for the 
field strength obtained in section I4~2| we have P mag > 1-96 x 10" 7 dyne cm -2 . This value is 
at least one order of magnitude higher than the radiation and ionization pressures . Even 
we add the thermal pressure to the calculations (Pther — 1-4 x 10~ 9 dyne cm" 



Vallee 



( 119871 )). the energy injected by the ionization front is still lower than the magnetic force. 



Therefore the expanding Hll region is not enough to compress the magnetic field lines into 
the observed geometrical configuration, and the bending is produced by the gravitational 
pulling. 



4.4. Multiple Outflows 



Previous works reported that FIR 5 has an associated (redshifted) unipolar and highly 
collimated outflow with a mass of ~ 4 M and density of ~ 10 2 cm -3 ( jSanders fc Willner 
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19851 ; iRicher et al.l Il992t IChandler fe Carlstroml Il996f ). However, a rather complex 



morphology was proposed by 



Cherninl (119961 ) as indicated by their interferometric (BIMA) 



and single dish (NRAO 12 m telescope) combined maps of CO (1 — > 0). In those, in 
addition to the unipolar lobe associated to FIR 5, there are two other redshifted features 
along the North-South direction and practically parallel to the FIR 5 molecular outflow, but 
neither of them associated with it. These two components were named nsl and ns2 and are 



detected at ulsr velocities between 15 and 25 km s 1 . ns2 is asso ciated wi t 
also identified a blueshifted feature, ewl, extending east of FIR 6 



h FIR 6. They 



Cherninl ( 119961 ) proposed 



that the brightest outflow component apparently powered by FIR 5A has a layered velocity 
structure. Their lower resolution combined molecular maps (~ 4"5) are dominated by an 
unipolar red lobe composed by two parallel outflows at lower velocities (~ 20 km s _1 in 
their Figure 1) which merge into only one at high velocities . This redder emission, referred 
as nsl in their paper, is more collimated than the lower velocity components and arises 10" 



west of FIR 5: A. 



or internal shocks (IChernin fc Masson 



he author suggests that the nsl outflow is widened by jet-wandering 



19951 ) and is powered by a deeply embedded and 



undetected source rather than FIR 5 A due to its misalignment with it. 

In this work, we offer a different interpretation for this scenario. The SMA CO 
(3 — > 2) m a ps ha ve an angular resolution a factor of 2 higher tha n the combined maps 



of 



Cherninl (119961 ). Contrary to the suggestion by 



Cherninl (119961 ). our maps show that 



this outflow is clearly powered by FI R 5A inst e ad of by a faint, undetected low-mass star. 

iin] 



The overall morphology described by 



Cherninl (119961 ) is also observed in our maps. The 



main difference is that we detect high velocity gas which is offset by 10" west of FIR 5 A, 
coinciding in position with the previously undetected FIR 5-sw dust condensation. Then, 
two interpretations can be derived from those features. Firstly, it is possible that all 
components are part of a single but velocity-layered outflow, the two low velocity lobes 
tracing the cavity where the highest velocity outflow is located. The presence of the high 
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velocity lobe not only at the center of the cavity but also displaced from it could suggest 
that the outflow is precessing. Alternatively, the presence of the FIR 5-sw dust source 
associated with the western red shifted lobe, and in particular at high outflow velocities, 
suggest that this lobe could be an independent molecular outflow. As in the case of FIR 
5A, this outflow would be also an unipolar outflow. 

Our CO (3 2) maps seem to indicate a possible interaction between the different 
outflows. In the blue lobe of Figure O there is extended emission centered in the dust 
condensation FIR 6n (using the LCGR02 nomenclature) with an EW orientation. The 
emission is associated with an unipolar outflow detected from ~ 1.0 to 7.8 km s _1 and is 
characterized by a wandering/wiggling morphology. The outflow suffers structural changes 
in its shape, represented by drastic rotations in PA. From being powered initially toward the 
NW direction, it assumes an almost horizontal distribution (PA = 94°) which corresponds 
to its brightest emission. Then, another structural change seen at 7.8 km s _1 results in an 
u— like shape. Both the bending from NW to the EW direction and the u— like structure 
coincide spatially with the projection of the red-lobe NS outflows powered by FIR 5A and 
FIR 5: sw. Figure [TT1 shows contours of velocity channels 15.6, 19.2 and 29.7 km s" 1 (black 
contours) overlaid with the 7.2 km s _1 channel (red contours) in the upper, middle and 
bottom panels, respectively. The variations in PA of the EW outflow may be somehow due 
to the interaction with the main outflow powered by FIR 5A and the high velocity emission 
from FIR 5-sw. Global inspection of these three panels tentatively leads to the hypothesis 
of a shock interaction between distinct outflows in such a way that the gas structure is 
modified. The spectrum exhibited in Figure H2] was obtained for the supposedly interacting 
zone of the panels of Figure [TTJ A box of -6" < Aa < -18" and -9" < AS < -24" was 
used to select a region where all three outflows components (EW, main and high velocity 
features) contribute to the emission. The lack of CO emission seen at ~ 25 km s -1 is 
probably due to the cavity previously mentioned and corresponds to the velocity interval 
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between low and high red velocity components. The blue emission is related to the EW 
outflow. 



4.5. Unipolar molecular outflow 



Previous observations of NGC 2024 fail to detect a blue c ounterpart for the bright 



outflow powered by FIR 5 A (IRicher et al. 



claim a bipolarity for this outflow ( panders fc Willner 



1992 



Chernin 



1985 



1 99 6[) . However, othe r studies 



Barnes et al. 



19891 ) . In all 



cases, the red shift lobe which is ejected toward the south has a brighter, more extended 
and collimated emission than the putative blue counterpart. These works all identify this 
feature as the main component of this outflow, making the morphology of the blue lobe, if 
it really exists, of unclear nature. 

This asymmetrical pattern of the outflow powe r ed by FIR 5A could be explained by 



the cloud morphology proposed by 



Matthews et al 



(120021 ) . illustrated in Figures 6 and 8 



in their paper. They indicate how the NGC 2024 Hn region could be seen from the west 
and north directions, respectively. In this scenario, the dust cores appear in the dense 
molecular cloud behind the Hn expanding front (using the line-of-sight as reference), with 
FIR 5 projected right below the interface zone. Consequently, the bipolar molecular outflow 
ejected from this core will have the blue molecular component destroyed by the UV photons 
produced in the Hn region, since it points right toward it, but the red lobe would remain 



intact. Despite the fact that 



Barnes et al. 



(119891 ) did not know accurately which core is the 



driving source of the supposedly blue outflow, these authors found that the total luminosity 
of the nebula is comparable to the flow energy. Therefore, some interaction between both 
could be expected. 
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5. Conclusions 

In this paper, we report SMA polarization observations of the intermediate-mass 
protostellar core NGC 2024 FIR 5. Data acquisition was done using the polarimetric 
capabilities of the SMA combined with wide spectral window receivers. The polarized flux 
appears distributed in three components: two of them around the peak of total intensity 
(Stokes I) and another component arising from the elongated portion of the core. The 
overall polarization portion resembles a partial hourglass morphology due to a possible 
ambipolar diffusion phenomenon taking place in the core. 

The magnetic field strength was estimated in 2.2 mG. The estimates of turbulent- 
to-magnetic energy and mass-to-flux ratio are consistent with a supercritical highly 
magnetized core. In previous works, magnetized collapsing cores were also observed in 
high-mass protostars. In general, ambipolar diffusion seems to affect core evolution globally, 
independent of the mass range. 

The absence of a symmetrical field morphology gives rise to different interpretations for 
the field structure in the core. The dust cores in NGC 2024 may be affected by an expanding 
ionization front compressing the molecular gas. It could be perturbing the field structure 
at smaller scales. The bended lines observed in our SMA maps could be the consequence 
of the radiation pressure of the hot component. Previous VLA 6 cm observations trace the 
foreground Hll region as an extended emission produced by the 02-B2 ionization source 
IRS 2b. However, our estimations of radiation pressure due to the expanding shell does not 
overcome the magnetic pressure generated by the field lines. So, the asymmetrical magnetic 
field is more likely due to depolarization effects arising in the position of the previously 
unresolved FIR 5B source. 

A complex outflow morphology was observed toward FIR 5. Several collimated features 
were detected toward FIR 5, FIR 6 and FIR 5-sw. We speculate about a possible flow 
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interaction between distinct components. It could explain the structural changes observed 
in some outflows. The brighter emission powered by FIR 5A has a clumpy structure and 
arises highly collimated in a NS orientation. The absence of a blue lobe counterpart can be 
attributed to the expanding Hll region to the north of the core. The UV radiation field 
may be responsible for dissociating the molecular structure of the outflow, destroying this 
component. 

FOA would like to thank all the staff at the SMA in Hawaii and Cambridge, MA. 
FOA and JMG are supported by the MICINN AYA2008-06189-C03 and the AGAUR 
2009SGR1172 grants. 
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Fig. 1.— VISTA image of the Flame Nebula (NGC 2024). The obscuring dust lane that 
exists foreground to the bright Hn emission is seen almost transparent in this near infra-red 
image. The glow of NGC 2023 and the Horsehead Nebula are seen in the lower portion of 
the image. 
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Fig. 2. — Dust continuum map of FIR 5 with quasi-uniform weighting (robust = — 1). 
Contours are drawn at —3, 3, 4, 6, 9, 13, 19, 25, 32, 42, 52, 62 a (1-er ~ 18 mJy beam -1 ). 
The half power beam width (HPBW) of the synthesized beam is l'/96 x l'/41 and the position 



angle is —70.6°. The crosses indicate the dust continuum sources detected by 
with BIMA. 
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Fig. 3. — Maps of Stokes Q {top panel) and U {bottom panel) emission. Dashed and solid thin 
contours correspond to negative and positive polarized emission, repectively. The contours 
start at — 2-cr and 2-cr level with steps of 1-er (1-er = 5.3 mJy beam -1 ). The absolute Q 
and U peak fluxes are 0.056 Jy beam _1 and 0.047 Jy beam -1 , respectively. The thick grey 
contours show the Stokes I emission. Crosses indicate the position of the two dust intensity 
peaks. The synthesized beam of the maps is shown in the bottom left corner of the bottom 
panel. 
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Fig. 4. — Contour map of the dust continuum emission overlapped with the linear polariza- 
tion vectors (black vectors) towards NGC 2024 FIR 5. Gray scale correspond to the polarized 
dust intensity. Contours levels are —3, 3, 4, 7, 11, 16, 22, 32, 42, 52, 62-x the rms noise of 
the dust emission (~ 19 mJy beam -1 ). The length of each segment is proportional to the 
degree of polarization (the length corresponding to a polarization level of 10% is indicated at 
the image bottom) . The synthesized beam is of 2'.'45 x l'/48 with a position angle of —40°. 
Vectors are sampled as 2/3 of a beam. 
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Fig. 5. — Histogram of position of polarization angles. The three polarized components of 
FIR 5 polarization map are clearly seen in this plot. 
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Fig. 6. — Channel maps of the CO (3 — > 2) emission associated to FIR 5 and FIR 6 dust cores. 
Contours levels are -4, 4, 8, 12, . . . to 36 x 0.25 mJy beam -1 (the rms noise of the map). The 
value of the Vlsr is shown in the top left corner of each panel. Source positions are indicated 
as crosses and labelled in the first panel. Magenta, green, red and brown arrows indicate 
the position of the FIR 6, FIR 5 A, the precessing and the FIR 5-sw outflows respectively. 
An ellipse indicates the supposedly cavity produced by the high- velocity components of the 
main FIR 5 outflow component. 
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Fig. 7. — Position- Velocity plot of the CO (3 — > 2) emission centered close to source FIR 5A 
(one arcsecond to the west) and along the North (positive offsets) to South (negative offsets) 
direction. Contours levels are —5, —3, 3, and then steps of 3 times 0.3 Jy beam -1 , the rms 
noise of the channel maps where the cut was obtained. The position of the driving source 
of the redshifted outflow, FIR 5A, is indicated with a dashed line. The spatial overlap with 
the East ouftflow associated with FIR 6 is also shown with a dashed line. 
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Fig. 8. — Distribution of polarization toward NGC 2024 FIR 5. Left panel: Polarization 
intensity versus total intensity. Right panel: Polarization intensity versus radius with respect 
to the peak of Stokes I emission. 
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Fig. 9. — Plane-of-sky field geometry for NGC 2024 FIR 5. Contours, beam and vector scale 
are the same than in Figure HJ Vectors are plotted at 2<jp level and relaxed to lop in the 
upper right corner. 
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Fig. 10. — VLA 6 cm emission from the Hn region in NGC 2024 (data from the VLA 
Archive). Grey scale-filled contours are 0.5, 1.5, 2, 3, 4, 5, 6, 7, 9, 11, 13.1a (la ~ 9 
mJy beam -1 ). The beam size of 8.6" x 7.5" and PA of 37° is shown in the lower left corner. 
Crosses indicate the location of FIR 5A and FIR 5B sources. The star indicates the position 
of the ionization source. 
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Fig. 11. — Superposition of contours from 15.6, 19.2 and 29.7 km s _1 velocity channels (black 
contours) over 7.2 km s _1 channel (red contours). Intensity contours are -3, 3, 4, 5, 7, 9, 12 
a (1 a ~ 0.57 mJy beam" 1 ). Source positions are maked as crosses. 
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Fig. 12. — Spectrum of the interacting zone between the EW outflow, powered by FIR 6, 
and the high velocity feature apparently powered by FIR 5-sw. The spectrum was obtained 
for a velocity range of to 33.2 km s _1 . The three peaks correspond to the emission from 
the EW outflow (the blue shifted peak at ~ 7 km s _1 ), the main lobe powered by FIR 5A 
at ~ 18 km s _1 and the high velocity lobe arising from FIR 5-sw at ~ 30 km s _1 . 



Table 1. Parameters of the continuum and line observations 



Observations Rest frequency HPBW PA a Spectral resolution Peak of emission rms noise 

(GHz) (arcsec) (°) km s -1 (Jy beam -1 ) (Jy beam -1 ) 



Continuum 345.8000 2.45 x 1.48 -39.8 1.19 0.019 b 

CO (3->2) 345.7960 2.87 x 1.67 -37.6 0.7 9.5 0.57 



a Position angles are measured from North to East. 

b The rms noise of Stokes I emission, obtained with a robust weight of 0.5 
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Table 2. FIR 5: main component 3 



Dust «2000 b 5 20 oo b Pcak b Total c FWHM Deconvolved 11 Deconvolved 

condensation of intensity flux Gaussian fit size PA 

(Jybeam -1 ) (Jy) (arcsec) (arcsec) (°) 

FIR 5A 05 41 44.258 -01 55 40.94 1.16(2) 2.44(6) 2.8 x 2.4 2.31(6) X 1.58(8) 47(5) 

FIR5B 05 41 44.510 -01 55 42.35 0.40(2) 1.26(8) 3.6 x 2.8 3.1(2) x2. 3(1) 130(11) 

Note. — Units of right ascension are hours, minutes and seconds, and units of declination are degrees, arcminutes, and 
arcscconds. 

a Fit error of the last digit in parenthesis. 

b Estimated using Miriad's "MAXFIT" task. 

c Values derived with AIPS's "IMFIT" task. 
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Table 3. Sub- millimeter dust condensations 



Dust 


^2000 b 


C2000 


Peak 


Flux 


condensation 






of intensity 13 


density 








(mJy beam -1 ) 


(mJy) 


FIR 5-sw 


05 41 43.667 


-01 55 49.05 


77(18) 


81(21) 


FIR 5-ne 


05 41 45.043 


-01 55 31.60 


131(18) 


202(30) 


FIR 6n a 


05 41 45.193 


-01 56 00.50 


70(18) 


40(15) 


FIR 6c a 


05 41 45.134 


-01 56 04.01 


124(18) 


116(22) 



a According to lLai et al.l ( 120021 ) numbering. 
Estimated using Miriad's "MAXFIT" task. 
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Table 4. SMA polarization data from NGC 2024 FIR 5 



A RA a A Dec P eP a P h 


Ip c 


d 


cO 


(arcsec) (arcsec) (%) (%) 


(Jy beam -1 ) 


(°) 


(°) 



5.4 


-2.1 


6.00 


2.0 


3.00 





,018 


35.723 


9.194 





-2.1 


13.4 


3.1 


4.32 





,027 


-34.364 


6.093 


5.4 


-1.5 


9.20 


2.8 


3.29 





,019 


27.954 


8.355 


4.5 


-1.5 


3.40 


1.4 


2.43 





.013 


24.122 


12.244 


1.8 


-1.5 


3.60 


1.6 


2.25 





.013 


-48.573 


12.838 


0.9 


-1.5 


3.30 


1.2 


2.75 





.016 


-47.342 


10.207 





-1.5 


5.20 


1.1 


4.73 





,027 


-35.171 


6.1 


-0.9 


-1.5 


8.00 


1.9 


4.21 





,025 


-34.346 


6.549 


1.8 


-0.9 


2.00 


1.0 


2.00 





,012 


-49.173 


13.955 


0.9 


-0.9 


1.70 


0.7 


2.43 





,014 


-55.723 


11.734 





-0.9 


1.80 


0.6 


3.00 





.017 


-39.46 


9.757 


-0.9 


-0.9 


3.40 


1.0 


3.40 





.019 


-37.513 


8.696 


1.8 


-0.3 


1.50 


0.7 


2.14 





.012 


-60.593 


13.363 


0.9 


-0.3 


2.10 


0.6 


3.50 





,022 


-64.888 


7.45 





-0.3 


1.40 


0.5 


2.80 





,016 


-52.705 


9.965 


-0.9 


-0.3 


1.60 


0.8 


2.00 





,012 


-47.955 


13.225 


3.6 


0.3 


4.80 


2.4 


2.00 





.012 


32.233 


13.901 


1.8 


0.3 


2.40 


0.7 


3.43 





.019 


-64.96 


8.479 


0.9 


0.3 


4.00 


0.5 


8.00 





.042 


-67.692 


3.88 





0.3 


2.80 


0.5 


5.60 





,031 


-58.781 


5.187 
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Table 4 — Continued 



A RA a 


A Dec 


P 


cP 




Ip c 


d 


cO 


(arcsec) 


(arcsec) 


(%) (%) 




(Jy beam -1 ) 


(°) 


(°) 


-0.9 


0.3 


1.70 


0.7 


2.43 


0.013 


-56.207 


12.44 


2.7 


0.9 


3.90 


2.1 


1.86 


0.011 


-10.944 


15.247 


1.8 


0.9 


3.90 


0.9 


4.33 


0.023 


-62.435 


6.924 


0.9 


0.9 


6.70 


0.6 


11.2 


0.061 


-70.852 


2.667 





0.9 


6.30 


0.7 


9.00 


0.051 


-63.767 


3.165 


-0.9 


0.9 


3.30 


1.0 


3.30 


0.018 


-58.229 


8.806 


1.8 


1.5 


6.80 


1.7 


4.00 


0.023 


-57.537 


7.051 


0.9 


1.5 


9.40 


0.9 


10.4 


0.061 


-73.461 


2.66 





1.5 


11.7 


1.1 


10.6 


0.063 


-69.462 


2.569 


-0.9 


1.5 


8.60 


2.0 


4.30 


0.025 


-61.584 


6.407 


0.9 


2.1 


11.1 


1.5 


7.40 


0.043 


-71.543 


3.729 





2.1 


15.0 


1.7 


8.82 


0.054 


-72.681 


2.981 


0.9 


2.7 


14.2 


3.1 


4.58 


0.028 


-61.117 


5.685 





2.7 


12.6 


2.7 


4.67 


0.029 


-71.458 


5.554 



a Offset respect to the peak of total intensity (same for dec- 
lination). 

b Signal-to-noise ratio of polarization. 
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c Polarized intensity xlO -2 . 

d Position angles are measured from North to East. 



